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ABSTRACT 



We present ID numerical simulations aimed at studying the hot-flasher scenario for the formation of He-rich subdwarf stars. 
Sequences were calculated for a wide range of metallicities and physical assumptions, such as the stellar mass at the moment of 
the helium core flash. This allows us to study the two previously proposed flavors of the hot-flasher scenario ("deep" and "shallow" 
mixing cases) and to identify a third transition type. Our sequences are calculated by solving simultaneously the mixing and burning 
equations within a diffusive convection picture, and in the context of standard mixing length theory. We are able to follow chem- 
ical evolution during deep-mixing events in which hydrogen is burned violently, and therefore able to present a homogeneous set 
of abundances for different metallicities and varieties of hot-flashers. We extend the scope of our work by analyzing the effects of 
non-standard assumptions, such as the effect of chemical gradients, extra-mixing at convective boundaries, possible reduction in con- 
vective velocities, or the interplay between difussion and mass loss. Particular emphasis is placed on the predicted surface properties 
of the models. 

We find that the hot-flasher scenario is a viable explanation for the formation and surface properties of He-sdO stars. Our results also 
show that, during the early He-core burning stage, element diffusion may produce the transformation of (post hot-flasher) He-rich 
atmospheres into He-deficient ones. If this is so, then we find that He-sdO stars should be the progenitors of some of the hottest sdB 
stars. 

Key words, stars: evolution — Stars: horizontal-branch — Stars: subdwarfs — Stars: mass-loss 



1. Introduction 

Hot subluminous stars are an important population of faint blue 
stars that can roughly be grouped into the cooler sdB stars, 
whose spectra display no or only weak helium lines, and the 
hotter sdO stars, which have higher helium abundances and can 
even be dominated by helium. Subluminous O, B stars were 
identified with stars populating the hot end of the horizontal 
branch (HE) of some globular and open clusters (de Boer et al. 
1997). While the sdB stars form a homogeneous spectroscopic 
class, a large variety of spectra is observed among sdO stars 
(Lemke et al. 1997). This diversity of spectra and the helium- 
enhanced surface abundances observed in many of these stars 
pose a challenge to our understanding of stellar evolution. As a 
consequence, some non-canonical evolutionary scenarios were 
proposed for their formation. Among them, the merger of two 
white dwarfs (Saio & Jeffery 2000) and a core helium flash af- 
ter departure from the red giant branch (i.e. "hot-flasher sce- 
nario") off'er the most promising explanations of their formation 
(Stroer et al. 2007, Moehler et al. 2007, and Hirsch et al. 2008). 
The second possibility is the subject of the present article. The 
hot-flasher scenario was proposed to explain the existence and 
characteristics of blue hook stars (such as helium or carbon sn- 
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hancement) of some globular clusters e.g. wCen and NGC 2808 
(D'Cruz et al. 2000, Moehler et al. 2002, Moehler et al. 2004, 
and Moehler et al. 2007). The coexistence of extreme horizon- 
tal branch (EHB) stars and helium core white dwarfs in globular 
clusters (Calamida et al. 2008) also agrees with a natural pre- 
diction of the hot-flasher scenario (see Castellani et al. 2006). 
Although in this scenario mass loss must be increased artificially, 
from values usually adopted in the modeling of low mass AGB 
stars, the driving mechanism behind RGB mass loss is not well 
understood (Espey & Crowley 2008). Differences in metallicity, 
initial He abundances or rotational velocities may, in principle, 
produce enhanced mass loss rates (the last two by an increase in 
the RGB-tip luminosity). The later is particularly interesting in 
view of the fact that some He-sdO stars are known to be rela- 
tively fast rotators (in comparison with sdB stars; Hirsch et al. 
2008). The existence of high mass-loss rates in old metal-rich 
populations is also an important ingredient of one of the most 
successful hypothesis for modeling the UV-upturn of elliptical 
galaxies (Yi & Yoon 2004). If this is so, it is very probable that 
a fraction of stars in these populations will undergo the helium 
core flash (HeCF) after departing from the first red giant branch 
(RGB) — i.e. will be "hot-flashers" as coined by D'Cruz et al. 
(1996). Finally, in high stellar density enviroments, it is possible 
that enhanced mass loss might result from star-to-star encoun- 
ters. 

In general, low mass stars undergo the HeCF at the tip of 
the RGB. However, Castellani & Castellani (1993) showed that, 
if sufficient mass is lost on the RGB, the star can depart from 
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the RGB and experience the HeCF while descending the hot 
white dwarf cooling track. Afterwards D'Cruz et al. (1996) ana- 
lyzed the frequency of these HeCFs at high temperatures ("Hot 
Flashers") by assuming different values of the mass-loss effi- 
ciency and suggested that this scenario might provide a path for 
populating the hot end of the HB. Sweigart (1997) noticed the 
similarity between these hot-fiashers and the late helium shell 
flashes that produce "bom-again AGB stars" (Iben 1984) and 
showed that — as happens in the born-again scenario — the con- 
vective zone, generated by the extremely high energy liberation 
rate (Lhb ~ 10 '^L©) during the flash, reaches the H-rich enve- 
lope and engulfs it. As a consequence, much or almost all H is 
burned leading to a He-enriched surface abundance. Brown et al. 
(2001) explored this scenario to explain anomalies in the hot HB 
of NGC 2808. However, it was not until the work of Cassisi et al. 

(2003) that simulations including the violent mixing and burning 
of H were performed and yields for the abundances of helium, 
carbon and nitrogen were calculated consistently. All of these 
works were focused mainly on globular clusters and therefore 
on low metalUcity progenitors. However, counterparts of these 
hot-flashers are also expected in the Galactic disk. This was con- 
sidered by Lanz et al. (2004) who assessed the possibility that 
(field) He-sdB stars could be formed by the hot-flasher scenario. 
In this work, the authors classified three different clases of hot- 
flashers:(l) early hot-flashers, which experience the HeCF dur- 
ing the evolution at constant luminosity in the HR diagram (af- 
ter departing from the RGB-tip) and become hot subdwarf stars 
with standard H/He envelopes; (2) late hot-flashers with shallow- 
mixing (SM), which become He-enriched hot subdwarfs due to 
convective dilution of the envelope (at the low-luminosity, low 
temperature region of the HR diagram) into deeper regions of the 
star; and (3) late hot flashers with deep-mixing (DM) in which 
the H-rich envelope is engulfed and burned in the convective 
zone generated by the primary HeCF. Due to the numerical diffi- 
culties related to the calculation of the violent H burning that oc- 
curs during a deep-mixing event (see Cassisi et al. 2003), Lanz 
et al. (2004) failed to provide surface abundances in this case. 
FinaUy, Stroer et al. (2005, 2007) and Hirsch et al. (2008) an- 
alyzed the possibiUty that field He-sdO could be generated by 
this mechanism and concluded that the idea could not be re- 
jected, although some serious discrepancies exist between mod- 
els and observations. However, the lack of a quantitative study of 
hot-flasher events and the chemical composition that arise from 
them for a wide range of metallicities, ensures that comparison 
with observations remains difficult. We note that only one pre- 
vious work (Cassisi et al. 2003) has calculated the abundances 
produced by the deep-mixing event but only for one remnant 
mass and metaUicity (Z=0.0015). On the other hand Lanz et al. 

(2004) calculated several sequences for solar-like metallicities 
but no calculation of the final surface abundances for the deep- 
mixing events were presented. The amount of observational data 
for these stars is increasing rapidly as a byproduct of surveys 
of faint blue objects such as SDSS and SPY. As concluded by 
Cassisi et al. (2003) to assess the viability of the hot-flasher sce- 
nario, a study of the dependence of the hot-flasher predictions 
(Teff, g and abundances) on the physical details is needed. 

The core feature of this article is to present a homogeneous 
set of simulations of the hot-flasher scenario for a wide range 
of cases in which the chemical abundances of the models are 
consistently calculated. In particular we present simulations for 
4 different metalUcities and different post-RGB remnant masses 
(which determines the luminosity of the model at the moment 
of the HeCF). The work is organized as follows. In the follow- 
ing section, we discuss the main ingredients of the simulations 



performed in this work. We then, in Sect. 3, discuss different 
flavors of hot-flasher and He-enrichment and describe our sim- 
ulations results under standard assumptions. In Sect. 4, we an- 
alyze deviations from the standard assumptions such as possi- 
ble extra-mixing at convective boundaries, the effect of chem- 
ical gradients, or a possible reduction in convective velocities. 
We also address the possible effects of element diffusion in the 
outer layers. In Sect. 5, we discuss the results of our standard 
sequences and compare them with observations. FinaUy, in Sect. 
6, we present some concluding remarks. 



2. Numerical and physical details 

The sequences presented in this work were calculated with 
LPCODE, a numerical code for solving the equations of stellar 
evolution, which was recently used to model the formation of 
H-deficient post-AGB stars in the born-again scenario (Miller 
Bertolami et al. 2006, and Miller Bertolami & Althaus 2007) 
and, with minor changes, to study global properties of the evolu- 
tion of low and intermediate stars (Serenelli & Fukugita 2007). 
LPCODE is a Henyey-type stellar evolution code designed specif- 
ically to compute the formation and evolution of white dwarf 
stars and described extensively by Althaus et al. (2003, 2005 and 
references therein). We therefore only briefly refer to the code to 
mention some particular features and some modifications of spe- 
cial interest to this work. 

Radiative opacities are those of OPAL (Iglesias & Rogers 
1996), including opacities for C-rich mixtures as expected for 
the regions with violent proton burning during a deep-mixing 
event (Cassisi et al. 2003) complemented at low temperatures 
with the molecular opacities of Alexander & Ferguson (1994). 
In the present work, conductive opacities for the degenerate He- 
core were included according to Cassisi et al. (2007). As a conse- 
quence, we note that the core masses of our models at the HeCF 
are sUghtly lower than in previous works with AMne-core ~ 
0.006 Mq for all metallicities (see Cassisi et al. 2007). Neutrino 
emission due to plasma processes was calculated as described 
by Haft et al. (1994). We used a nuclear network that considers 
16 elements and 34 nuclear reactions of the p-p chains, CNO 
bi-cycle, helium burning, and carbon ignition that are identical 
to those adopted by Althaus et al. (2005) with the exception of 
^^C+p ^i^NH-y -^^^C+e+ + v, and ^^C(p,yy'^N, of special 
interest in this work, which are now taken from Angulo et al. 
(1999). 

In this work, mixing and burning of elements are solved si- 
multaneously in the context of diffusive convective mixing. This 
is performed by solving the set of equations 



dX 
dt 



dm 



/ , \2 dX 



(1) 



where Z is a vector containing the abundances of all consid- 
ered elements (see Althaus et al. 2003 for details of the nu- 
merical procedure). The efficiency of convection (or any other 
mixing process) is described by adopting the appropriate diffu- 
sion coefficient, D. In the present work convection was solved 
(generally, see Sect. 4) with the standard mixing length theory 
(MLT), adopting a value of a = 1.61 for the free parameter of 
the MLT. This value aUowed us to reproduce the present lumi- 
nosity and effective temperature of the sun, logTeff = 3.7614 
and Lq = 3.842 x lO^^erg s~\ at an age of t = 4570 Myr, when 
adopting Z - 0.0164 and X = 0.714 in agreement with the Z/X 
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Mass at ZAMS 
(Mo) 



Initial Abundances 
(X/Y/Z) 



0.88 0.769/0.230/0.001 

0.98 0.736/0.254/0.010 

1.03 0.702/0.278/0.020 

1.04 0.668/0.302/0.030 

Table 1. Initial (ZAMS) values of our standard sequences. 



value of Grevesse & Sauval (1998). In the context of the MLT, 
convection is described by 



D 



^(1 

Kp 



•yS)V,d(V„d-V) 



1/3 



(2) 



Initial (ZAMS) masses for our simulations (Table 1) were 
chosen to obtain ZAHB ages close to those of globular clusters 
(~ 12.5 Gyr, Salaris & Weiss 2002). However we do not expect 
initial mass to play a role in these simulations and the results 
should therefore be valid for sequences of ages at the moment 
of the HeCF similar to that of the disk (~ 8 Gyrfl Mass loss in 
LPCODE is now included with the prescription of Schroder and 
Cuntz (2005). However, as the details of mass loss in the RGB do 
not coiTespond to distinctive features in the models after the star 
has departed from the RGB, we enhanced mass loss artificially 
close to the RGB tip to obtain different H-rich envelope masses 
at the moment of the HeCF. Mass loss was stopped arbitrarily 
when the star reached Teff ~ 12 600 K. 

The initial helium content of the sequences (which span a 
wide metallicity range) was chosen to be 



Y = 0.23 + 2.4 X Z, 



(3) 



allowing comparison with previous works, that have adopted 
Y - 0.23 for metal-poor globular clusters, while at the same 
time is consistent with our adopted solar values and similar to 
present determinations of the chemical evolution of the galaxy 
(Flynn 2004, Casagrande et al. 2007). 



3. Description of standard sequences 

We performed simulations of the late hot-flasher scenario in 
which the surface composition is altered by both dilution and 
burning of the H-rich envelope. Evolutionary sequences were 
followed from the ZAMS, through to the hot-flasher event, and 
to the post horizontal branch evolution. Our standard sequences 
include neither overshooting at any convective boundary nor any 
other extra-mixing processes (such as semiconvection). In these 
sequences, we also neglected the effect of the /v-gradients at the 
H-He transition and convective velocities were those predicted 
by the MLT. Since our aim is to provide a complete and homo- 
geneous set of observable properties of the predictions of these 
events, we performed full evolutionary simulations for the differ- 
ent flavors of the hot-flasher scenario and four different metallic- 
ities (Z=0.001, 0.01, 0.02 and 0.03), and for He-core flashes at 
different points of the post-RGB evolution, i.e. at different rem- 
nant masses. Initial abundances and masses of these sequences 
are shown in Table 1 . 

Before describing of the surface properties of these se- 
quences, we provide a brief description of how H-deficiency is 
achieved in the simulations. 



' For example at Z=0.02 the mass of the He-core at the HeCF is less 
than 0.5% smaller for a higher mass progenitor (1.25 M©) that under- 
goes the HeCF ~ 6 Gyr after the ZAMS. 
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Fig. 1. HR diagrams for the three kinds of late hot-flashers de- 
scribed in this work. Sequences coiTespond to Z=0.02. For com- 
parison, an early hot-flasher sequence of 0.47426 Mq that fin- 
ishes with normal H abundances (broken lines), and a HeCF- 
manque sequence of 0.4607 Mq that ends as a He core white 
dwarf (dotted line) are also plotted, [color figure only available 
in the electronic version] 



Late hot-flashers, i.e. those that experience He-enrichment of 
their stellar surfaces, were classified by Lanz et al. (2004) into 
"shallow" and "deep" mixing cases. In shallow-mixing (SM) 
cases, no H-burning occurs and He-enrichment is attained only 
by dilution of the remaining H-rich envelope into deeper regions. 
In contrast, in deep-mixing (DM) cases the H-rich envelope is 
engulfed by He-flash convective zones and H is burned violently 
in the hot interior of the star. This picture is not entirely consis- 
tent with our results in particular because there are intermediate 
cases in which, although He-enrichment is due to dilution some 
inner parts of the H-rich envelope are burnt during the episode. 
In these cases L//, can be momentaneously as high as 10^ L© 
but only a fraction of the H-content is finally burnt. Since there 
is no sharp separation between these intermediate cases and the 
"shallow mixing" cases of Lanz et al. (2004), we retain the "shal- 
low mixing" denomination for both cases. However, we denote 
the existence of some important H-burning during the flash by 
labeling those cases in which H-burning luminosity reached val- 
ues momentaneously higher than 10"* Lq as SM* (see Table 2). 
On the other hand, the separation between DM and SM cases is 
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clearly characterized by the development of a violent runaway 
burning of protons during DM cases and, especially, by the al- 
most complete burning of the original H content of the remnant. 
In the following paragraphs, we briefly describe these different 
cases of late hot-flashers. 
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Fig. 2. Kippenhahn diagram of a deep-mixing case. The case 
plotted corresponds to a 0.47017 Mq sequence of initially 
Z=0.02. Red (black) points mark the upper (lower) convective 
boundaries of each model. Main convective zones have been 
shaded for clarity. Blue (green) zones mark the regions in which 
energy generation per unit mass is above e > 10^ erg g"' s"' 
(e > 10^ erg g ' s ' ). Dotted line marks the location of the inner 
boundary of the H-rich envelope, defined hy X ~ 0.001. [color 
figure only available in the electronic version] 



Since the evolution of convective zones during all cases of 
late hot-flashers have some common features, we discuss these 
first. First, as a consequence of neutrino cooling, the primary 
HeCF develops off-center A growing convective zone develops 
in the He-core above the point of maximum release of energy, 
with its upper boundary becoming progressively closer to the 
H-He transition. When the flash is fully developed and the struc- 
ture of the star changes as a consequence of the sudden energy 
released by the primary HeCF — the star is moving to low lumi- 
nosity and low temperature regions of the HR diagram see Fig. 
[1] — a very shallow inwardly growing convective zone develops 
due to opacity and lower surface temperatures of the models. It 
is this convective zone that ultimately exposes the H-deficient 
material when it merges with some of the convective zones de- 
veloped in the interior of the star. This is true for all kinds of late 
hot-flashers (see Figs.|2l[3] andHJi. 

The evolution of the DM episode is shown in Figs. [T] and |2] 
First, the He-flash starts off-center and similar to those at the tip 
of the RGB, developing a convective zone whose upper bound- 
ary becomes closer to the H-He transition as the He-burning 
luminosity rises. As reported by Sweigart (1997), the entropy 
barrier (due to previous CNO burning) is insufficiently high — 
because the H-burning shell is almost extinguished at the mo- 
ment of the flash — and the HeCF driven convective zone is able 
to penetrate the H-rich envelope, bringing protons into the C- 
enriched (a few percent by mass) He-burning convective zone. 
As a consequence, H burns producing an extra energy release. At 
some point, the rate of H-rich material brought down and burnt in 
the interior is so high that the energy released by H-burning be- 
comes very important and two important events occur. First, due 



to the rise in the entropy of the material at the location at which 
protons burn these layers become convectively stable splitting 
the previous convective zone into two, one powered by the He- 
flash and the other by H-burning. Second, the energy released 
by proton burning becomes the main source in driving further 
expansion the convective region leading to an unstable situation: 
a further increase in the convective zone produces an increase 
in the power released by H-burning, which leads to a further 
outward excursion of the convective boundary. This produces a 
H-flash in which almost all H-content of the remnant is rapidly 
burnt, an episode that is also characterized by a rapid outward ex- 
cursion of the outer boundary of the convective zone. This addi- 
tional growth in the convective zone can be clearly seen in Fig.|2] 
about 1 year after the primary HeCF. During the episode, the star 
moves to the low luminosity, low temperature region of the HR 
diagram (see Fig.[Tli, where a very shallow outer convective zone 
develops. This convective zone becomes progressively deeper 
and ultimately, some years after the primary HeCF, merges with 
the inner convective zones and the surface of the star becomes 
H-deficient (Fig.|2]i. After this, the star proceeds through several 
secondary helium flashes, like in the case of "canonical helium 
core flashes" (i.e. at the RGB tip), before finally settling at the 
ZAHB as an extremely H-deficient star (Xh ~ 10 - 10 *"). The 
hydrogen content decreases from a few IO'^Mq to usually less 
than 10'^ Mq (see Table 2) after the primary helium flash. 
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Fig. 3. Kippenhahn diagram of shallow-mixing with slight CNO 
burning (SM*). References are similar to those of Fig.|2] [color 
figure only available in the electronic version] 



If the flash occurs when the H-burning shell is somewhat 
more luminous (a SM* case), then the convective zone may 
come into contact with the inner part (with low H content) of 
the H-rich envelope, but the material that is brought downwards 
is never sufficient for the H-flash to develop. This is shown in 
Fig. |3] (see also the middle panel of Fig. [T]!. Then only a tiny 
fraction of the envelope is burned — e.g. the total H-content of 
the star changes, after the entire series of HeCFs, from 2.5 x 10 '* 
to 6.9 X 10"^ Mo, in the case shown in Fig.|3] and in the Z=0.001 
SM* case of Table 2 only about ~ lO"'* Mq is burnt during 
the primary HeCF. Afterwards, the outer boundary of the flash 
driven convective zone penetrates slowly, diluting the inner re- 
gions of the H-rich envelope until it finally merges with the in- 
wardly growing surface convective zone, turning the star into a 
hydrogen deficient star with low H content in the envelope {X of 
10"^ to 10""^ by mass fraction, see Table 3). 
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Fig. 4. Kippenhahn diagram of shallow-mixing without slight 
CNO burning. References are similar to those of Fig. |2l [color 
figure only available in the electronic version] 
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Fig. 5. Mass and metallicity of the standard sequences calculated 
for this work. The orange zone indicates the sequences that un- 
dergo a change in their surface abundances as a consequence 
of mixing and/or burning of the H-rich envelope. Dashed line 
marks the separation between DM and SM cases, [color figure 
only available in the electronic version] 



Finally, if the flash occurs when the H-burning shell is even 
more active (see Fig. |4] and also the top panel of Fig. [TJ, then 
the convective zone does not come into contact with the H-rich 
envelope, due to the entropy, barrier and no H is burned (this 
is the case originally referred to as "shallow mixing" by Lanz 
et al. 2004). After the He-flash, a small tail of the HeCF con- 
vective zone remains for many thousand years moving outwards 
in mass. The convective zone reaches the H-rich envelope and 
dilutes its inner parts into the H-free interior. Finally, it merges 
with the inwardly growing surface convective zone, producing a 
mild H-deficiency in the stellar surface. 

Results of numerical simulations of all these events are pre- 
sented in Table 2, where the timescales and H-content at different 
stages are given. 

3.2. Sequences in the mass-metallicity plane 

In Fig. |5] we show the standard sequences calculated for this 
work in the mass-metallicit}|3 plane. As previously mentioned 
the masses of hot-flasher remnants are slightly lower than in 
previous works due to the use of the new conductive opacities 
presented by Cassisi et al. (2007). 

Our results support the claim by Lanz et al. (2004) that 
shallow-mixing is more important in sequences of high metal- 
licity (see Fig.lU. While shallow-mixing accounts for only ~4% 
of the remnant mass range of late hot-flashers in lower metal- 
licity sequences (Z=0.001 and Z=0.01), it corresponds to 14- 
17% of the remnant's mass range in the Z=0.02 and Z=0.03 se- 
quences. The remnant mass range for hot-flashers is found to be 
slightly dependent on metallicity, varying from AM ~ O.O14M0 
for Z=0.03 to AM ~ O.OIOSMq for Z=0.001, although this in- 
crease is mainly due to the increase in the SM mass range, while 
the DM mass range is almost constant. It is not easy to trans- 
late this result into a mass range for initial masses because this 
depends strongly on the mass-loss prescription or mechanism. 



^ By metallicity we mean the original metal (ZAMS) content of the 
star. During this work we choose to label the sequences with their final 
mass instead of the mass loss rate — as done in previous works — be- 
cause it is the final mass of the remnant, and not the particular way in 
which mass is lost, which plays an important role in the development of 
the different kinds of hot-flashers. 



3.3. Surface properties 

In Table 3, we present the surface abundances of our simu- 
lations for selected sequences. We first note that, while DM 
episodes lead to surface hydrogen abundances preferentially be- 
tween 10"^ and 10"^ by mass fraction (although they can be as 
high as a few lO^'* in some simulations) the SM/SM* episodes 
offer a wide range of surface H abundances ranging from al- 
most solar to a few 10""' by mass fraction. The SM/SM* sce- 
narios therefore provide the intermediate abundances needed 
by Unglaub (2005) to reproduce He-sdB/He-sdO surface abun- 
dances in the light of calculations that include both the effects 
of diffusion and weak mass loss rates (M < lO^'^Mo/yr). This 
result is quite interesting because some effects not included in 
these simulations (particularly the effect of the Vyu-barrier) may 
increase the frequency of SM to DM episodes due to the intro- 
duction of an extra barrier against the penetration of the HeCF 
convective zone into the H-rich envelope (see Sect. 4). 

For the other element abundances, some trends are evident 
for the DM cases (see Table 3). In particular, '^C decreases by 
a factor of 2 from earlier (i.e. higher mass remnants) to later 
cases (i.e. lower mass remnantsfl This is because, for later DM 
episodes, the H-flash and the subsequent splitting of the convec- 
tive zone occur more rapidly than in earlier DM cases (see ' At 
values in Table 2) and less ^^C is therefore created by the time of 
the H-flash. When the convective zones split, the increase in '^C 
inside the (inner) He-flash driven convective zone does not affect 
the final surface abundances. On the other hand, '^C shows no 
definite trend and is almost constant for a given metallicity, and 
decreases slightly when metallicity is increased. In Table 3, we 
can see that ''*N clearly increases for later deep-mixing episodes 
(for a given metallicity). This can be understood easily as a con- 
sequence of the earlier development, in those cases, of the H- 
flash. Since there is a lower amount of in the convective 
zone when both splitting occurs and H-burning develops, and 
the H mass engulfed is almost constant (see Table 2), then there 
are more protons per each '^C nuclei and more '"^N is therefore 

^ This definition of "later" and "earlier" does not refer to the time 
between ZAMS and HeCF but to the location in the HR-diagram in 
which the flash takes place. 
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Final Mass 


Zo 


maximum 


maximum 


'Af 






'M« 


'M„ 




[Mo] 




Log(LHe/Lo) 


Log(LH/Lo) 


[yr] 


[Myr] 


[Myr] 


[Mo] 


[Mo] 


[Mo] 


0.49145 (SM*) 


0.001 


10.23 


7.44 


9006 


1.53 


66.73 


4.19x10"'' 


3.21x10"" 


3.21x10"" 


0.49104 (DM) 


0.001 


10.22 


9.24 


106.78 


1.55 


67.33 


4.17x10"'' 


2.92x10"^ 


2.29x10"* 


0.48545 (DM) 


0.001 


10.10 


10.09 


7.35 


1.55 


65.20 


4.30x10"'* 


2.37x10"' 


4.16x10"'* 


0.48150 (DM) 


0.001 


10.09 


10.71 


4.11 


1.57 


69.60 


4.41X10-* 


2.10x10"'^ 


2.46x10"' 


0.47770 (SM) 


0.01 


9.83 


3.58 


79123 


1.71 


71.82 


2.93x10"" 


2.44x10"" 


2.44x10"" 


0.47744 (SM*) 


0.01 


9.83 


6.91 


18471 


1.66 


71.96 


2.93x10"" 


1.38x10"" 


1.31x10"" 


0.47725 (DM) 


0.01 


9.82 


9.12 


2960 


1.71 


72.69 


2.92x10"" 


2.33x10"* 


1.85x10"' 


0.46921 (DM) 


0.01 


9.66 


9.91 


54.7 


1.71 


84.69 


3.06x10"" 


6.72x10"' 


8.35x10"'* 


0.46644 (DM) 


0.01 


9.71 


9.70 


44.1 


1.65 


78.76 


3.11x10"" 


1.96x10"' 


3.78x10"'* 


0.47378 (SM) 


0.02 


9.71 






1.80 


73.89 


2.62x10"" 


2.62x10"" 


2.61x10"" 


0.47250 (SM*) 


0.02 


9.69 


5.78 


3053 


1.80 


78.99 


2.54x10"" 


9.66x10"^ 


6.85x10"^ 


0.47112 (DM) 


0.02 


9.66 


9.33 


706 


1.82 


76.48 


2.55x10"" 


8.27x10"' 


8.95x10"'* 


0.46410 (DM) 


0.02 


9.52 


10.06 


135.8 


1.83 


88.10 


2.66x10"" 


7.32x10"' 


7.91x10"'* 


0.46150 (DM) 


0.02 


9.60 


10.25 


132.4 


1.74 


81.75 


2.70x10"" 


2.01x10"' 


3.48x10"'* 


0.46521 (SM) 


0.03 


9.45 






1.87 


83.33 


2.67x10"" 


2.67x10"" 


2.56x10"" 


0.46470 (SM) 


0.03 


9.45 


2.67 


78688 


1.89 


79.70 


2.48x10"" 


1.70x10"" 


3.44x10"'^ 


0.46367 (SM*) 


0.03 


9.43 


5.80 


17248 


1.90 


80.76 


2.49x10"" 


7.23x10"^ 


1.16x10"' 


0.46282 (DM) 


0.03 


9.41 


9.00 


10926 


1.93 


81.47 


2.49x10"" 


4.41x10"* 


3.28x10"' 


0.45660 (DM) 


0.03 


9.27 


9.80 


975 


2.00 


87.76 


2.49x10"" 


7.05x10"' 


5.52x10"'* 


0.45234 (DM) 


0.03 


9.31 


9.49 


1055 


1.94 


93.93 


2.58x10"" 


5.64x10"' 


6.47x10"'* 



Table 2. Properties of some of the standard sequences calculated for this work. Fifth, sixth, and seventh column indicate the interval 
of time between the maximum of Log(LHe/i0) and hogiL-nj Lq) (' Af), from the primary He flash to the ZAHB (^Af) and from the 
ZAHB to the end of the EHB (TAHB;^Af). The final three columns provide the total H content of the sequence just before the He 
flash ('Mh), after the primary He core flash CMh) and at the ZAHB (^M^). 



Final Mass [Mo] 


Zo 


H 


He 


'^C 




N 


O 


0.49145 (SM*) 


0.001 


0.2137 


0.7475 


0.0365 


1.57 X 10"* 


1.22 X 10"" 


1.52 X 10- 


-3 


0.49104 (DM) 


0.001 


2.42 X 10"" 


0.9450 


0.0267 


7.47 X 10"^ 


0.0106 


4.77 X 10" 


-5 


0.48545 (DM) 


0.001 


4.60 X 10"* 


0.9524 


0.0264 


7.73 X 10"3 


0.0134 


4.53 X 10- 


-5 


0.48150 (DM) 


0.001 


2.41 X 10"* 


0.9666 


0.0107 


5.79 X 10"3 


0.0185 


3.79 X 10- 


-5 


0.47770 (SM) 


0.01 


0.0780 


0.8682 


0.0426 


6.26 X 10"* 


1.30 X 10"' 


6.84 X 10- 


-4 


0.47744 (SM*) 


0.01 


0.0261 


0.8996 


0.0642 


2.36 X 10"* 


6.59 X 10"" 


7.43 X 10" 


-4 


0.47725 (DM) 


0.01 


2.27 X 10"' 


0.9464 


0.0382 


6.12 X 10"^ 


4.35 X 10"' 


2.94 X 10- 


-4 


0.46921 (DM) 


0.01 


1.02 X 10"' 


0.9560 


0.0211 


6.13 X 10"^ 


0.0128 


2.92 X 10- 


4 


0.46644 (DM) 


0.01 


4.61 X 10"* 


0.9592 


0.0193 


5.52 X 10"3 


0.0122 


3.16 X 10- 


-4 


0.47378 (SM) 


0.02 


0.3822 


0.5978 


1.23 X 10"' 


8.96 X 10"' 


7.57 X 10"' 


4.69 X 10" 


-3 


0.47250 (SM*) 


0.02 


9.56 X 10"-'' 


0.9273 


0.0425 


1.93 X 10"' 


2.94 X 10"' 


8.64 X 10" 


-4 


0.47112 (DM) 


0.02 


9.05 X 10"* 


0.9389 


0.0420 


4.27 X 10"^ 


4.88 X 10"' 


6.96 X 10- 


-4 


0.46410 (DM) 


0.02 


1.05 X 10"' 


0.9460 


0.0273 


6.41 X 10"^ 


0.0120 


7.680 X 10 


-4 


0.46150 (DM) 


0.02 


4.66 X 10"* 


0.9480 


0.0263 


5.67 X 10"3 


0.0116 


1.01 X 10" 


-3 


0.46521 (SM) 


0.03 


0.2227 


0.7271 


0.0197 


7.79 X 10"' 


8.09 X 10"' 


5.18 X 10" 


-3 


0.46470 (SM) 


0.03 


4.71 X 10"^ 


0.9239 


0.0380 


4.12 X 10"" 


6.10 X 10"' 


1.58 X 10" 


-3 


0.46367 (SM*) 


0.03 


1.63 X 10"3 


0.9249 


0.0390 


3.25 X 10"3 


6.48 X 10"' 


1.62 X 10 


-3 


0.46282 (DM) 


0.03 


7.48 X 10"* 


0.9333 


0.0420 


1.74 X 10"3 


7.04 X 10"' 


1.55 X 10 


-3 


0.45660 (DM) 


0.03 


9.22 X 10"* 


0.9382 


0.0341 


4.09 X 10"^ 


0.0103 


1.67 X 10 


-3 


0.45234 (DM) 


0.03 


1.94 X 10"' 


0.9411 


0.0273 


4.76 X 10"3 


0.0134 


2.14 X 10- 


-3 



Table 3. Surface abundances (by mass fractions) of some of the sequences calculated in this work. 



formed. Then, at a given metallicity, the ratio of C/N can vary 
by more than a factor of two, depending on the moment of the 
post-RGB evolution at which the deep mixing occurs. This is 
interesting in view of the two classes of He-sdO stars (N- and 
C- rich) presented by Stroer et al. (2007). The dependence of 
the C/N ratio as a function of initial metallicity indicates that 
the lower the metallicity, the higher the N abundance. However, 
no clear correlation appears to be present between the location 
of our sequences with higher N/C ratios and the N-rich class of 
Stroer et al. (2007) (see Fig. IS. 

In the case of shallow-mixing events, '^C is orders of mag- 
nitude lower than for DM because no significant (if any) amount 
of H burns during the flash. At first, it may appear strange that N 



abundances in some cases, do not differ significantly (although 
lower) from those of DM episodes. However, during SM events, 
the material dredged up to the surface corresponds to material 
previously processed by the CNO cycle in the RGB, which trans- 
forms almost all primordial C and O into N. In this case, the N 
present in the surface is therefore strongly dependent on the ini- 
tial metallicity of the stars. This explains why the N abundance 
of SM events in Table 3 show a strong dependence on Zq. The 
final ''C abundances in SM depends mainly on the interplay be- 
tween the enormous amount of '^C created by the primary He 
core flash and the intensity of the dilution process when the in- 
ner and outer convective zones merge. 
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As noted by Stroer et al. (2007), although the tracks shown in 
Fig.|6]pass close to the location of He-sdO stars derived by Stroer 
et al. (2005, 2007), their distribution in the g-T^ff diagram does 
not coiTespond to the one expected from the simulations. From 
the point of view of the models, stars should cluster towards the 
location of the EHB, 4.5 < log^eff < 4.65 and 5.9 < logg < 6.2, 
because the evolution from the primary HeCF to the settlement 
in the ZAHB is far more rapid (occuring in less than 2x10^ yr) 
than the He-core burning stage (65 - 90 x 10^ yr, see Table 2). 
This clustering of He-sdO stars is not observed in the infeiTed 
g-Teff values of real stars (see Fig. |6]l. A possible explanation 
for this discrepancy may be that helium-rich atmospheres be- 
come hydrogen-rich, due to diffusion processes, in a timescale 
shorter than the duration of the He-core burning stage (EHB). In 
Sect. 4.4, we show that such a conversion is, in fact, possible. 
According to this scenario, late hot-flashers have a He-rich sur- 
face composition only during the evolution towards the EHB and 
then become stars with hydrogen dominated atmospheres. This 
imphes that He-sdO stars evolve into hydrogen-rich sdB/sdO 
stars, so that no clustering of He-sdO stars around the EHB is 
expected. This scenario is discussed in Sects. 4.4 and 5. 

When comparing theoretical predictions with the results of 
Stroer et al. (2007), we recall that surface properties in that work 
were obtained by adopting H/He atmospheres, which was found 
by Lanz et al. (2004) (who adopted H/He/C atmospheres in their 
analyzes) to be a possible cause of systematic effects. The ef- 
fective temperature inferred in real stars may also be systemati- 
cally overestimated due to line-blanketing effects (see Stroer et 
al. 2007). 

Finally, we note that the surface characteristics of our solar 
metallicity DM sequences agree closely with the values infeiTed 
by Lanz et al. (2004) for PG1544+488. This result clearly con- 
firms the statement of Lanz et al. (2004) that this object could 
be the result of a DM event. Our results, also, confirm that the 
surface abundances derived by these authors for JL87 agree with 
those predicted by the SM kind of hot-flasher. In contrast, no 
sequence presented abundances close to those of LB 1766, con- 
firming that this object may have had a very different history. 

3.4. Comparison with previous works 

We compare our results with the very few consistent simulations 
available in the literature. In the DM, we compare the 0.49125 
sequence of Cassisi et al. (2003), for which the He-core flash 
ocurrs at logL/L© ~ 0.79 (a relatively "early" deep-mixing 
episode in the terminology of Sect. 3.2) with the most similar 
models in our grid. This is shown in Table |4] where the results 





this work 


this work 


Cassisi et al. 


Mass [Mq] 


0.49104 


0.48545 


0.49125 


Zo 


0.001 


0.001 


0.0015 


'A/ 


107 d 


7.35 d 


150 d 


max. LhJLq 


1.66x10'° 


1.26x10'" 


3.7 xlO'" 


max. Lh/Lq 


1.74 xlO' 


1.23 xlO'° 


8.2x10' 


Mh burnt [Mq] 


3.88 xlO-"* 


4.30x10-'' 


3.7 xlO-" 


surf. H 


24x10-" 


4.6x10-'^ 


4 xlO-'' 


surf. He 


0.945 


0.952 


0.96 


surf. C 


0.0342 


0.0341 


0.029 


surf. N 


0.0106 


0.0134 


0.007 


surf. O 


4.8 xlO^^ 


4.5 xlO-"^ 


3.5 xlO-"^ 



Table 4. Comparison of characteristics of DM sequences of this 
work with the DM sequence of Cassisi et al. (2003). 



for the 0.48545 Mq and 0.49104 Mq (with Z=0.001) sequences 
of this work are compared with the DM sequence of Cassisi et 
al. (2003). We note in particular the good agreement between 
the earlier deep-mixing sequence (0.49104 Mq) and Cassisi et 
al. (2003) values, regardless of the different metallicity and mass 
in both sequences. 

For the SM case, the only available sequence is one of 
solar metallicity discussed by Lanz et al. (2004), which we 
can compare with our 0.47378 Mq, Z=0.02 SM sequence. In 
their sequence the shallow-mixing event (the merging of the 
inner and outer convective zones) occurs about 24000 yr af- 
ter the primary He-core flash compared with the 14000 yr in 
our sequence (see Fig. 3). The surface abundances in their se- 
quence are He/C/N=0. 5/0. 008/0. 005 qualitatively similar to ours 
(He/C/N=0.6/0.001/0.008). We can therefore conclude that there 
is qualitative agreement between Lanz et al. (2004) and our re- 
sults. 

4. Alternative sequences 

In this section we discuss the effect of some deviations from the 
standard assumptions made in the previous section. In particu- 
lar, we address the effect of possible extra-mixing at convective 
boundaries, different mixing velocities, and chemical gradients. 
Finally, we address the effect of diffusion in the outer layers of 
the standard sequences. 

4. 1 . Effect of tlie location of proton burning (aitering D) 

The calculation of ID-stellar evolution models through the H- 
flash in previous simulations (Sect. 3) were performed within the 
picture of diffusive convective mixing. This approach was used 
in the past for similar situations like late helium-shell flashes 
(i.e. born again episodes) experienced by post-AGB stars at the 
beginning of their white dwarf stage (e.g. Iben & Mac Donald 
1995). This picture, which implies that material is homoge- 
neously mixed in each mass shell, is very probably far from re- 
ality (see Woodward et al. 2008 for a discussion of this issue) in 
which convection appears to be an advection process. We now 
analyze how deep-mixing episodes can be altered if the loca- 
tion of maximum energy generation is displaced from the pre- 
dictions of MLT plus diffusive mixing. This is done by altering 
the value of the convective diffusion coefficient D in the mixing 
equations. We note, however, that because we still assume hor- 
izontal homogeneity, these results can only be considered poor 
approximations of the problerrQ- 

Surface abundances in simulations with reduced mixing ef- 
ficiencies are shown in Table 5. Our results are qualitatively 
consistent with the preliminary results quoted by Cassisi et al. 
(2003). They indicate, however, that a reduction in D of only 
two orders of magnitude is sufficient to produce an increase in 
the residual amount of H after a DM episode of about an or- 
der of magnitude. In view of the results of Unglaub (2005), it is 
interesting to note that a reduction in D of a factor of between 
1000 and 10000, produces H-abundances after the DM episode 
that would finally convert He-sdO stars into He-deficient sdB 
stars by gravitational settling on the EHB, even in the presence 
of weak winds (M ~ lO^'^Mg/yr). If major proton burning is 
shifted deeper into the star (Dx 10) then the extra energy released 
by proton burning is locally less relevant — in comparison with 

"* Horizontal inhomogeneities in '^C, '^C and H would certainly af- 
fect the final '''N yield. This is exactly the situation expected during the 
hydrogen flash (see Woodward et al. 2008). 
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0.47725 M^^ 


„ (DM) 


0.46644 M_^ 
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▲ 
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T 





Z=0.01 



He-sdO (C-rich) 
He-sdO (CN-class) 
He-sdO (N-iich) 
He-sdB 
sdB 

normal H composition 
0.46521 M^^^ (SM) 
0.46367 M^^^ (SM») 
0.46282 M_^^ (SM) 
0.45234 M (DM) 




Z=0.03 

.1.1. 



4 

4.5 
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6 

6.5 
4 

4.5 

5 ^ 
5.5 J 

6 

6.5 



5.1 5 4.9 4.8 4.7 4.6 4.5 4.4 



4.9 4.8 4.7 4.6 4.5 4.4 4.3 



Fig. 6. Values of T^fi and g of our sequences for different metallicities. Values inferred for He-sdO and sdB stars are from by Stroer 
et al. (2007) and Lisker et al. (2005), respectively. He-sdB data is taken from JefFery (2008) [color figure only available in the 
electronic version] 



Diff. Coef. 


H 


He 




'3C 


N 


o 




3M„ 


'Dmlt X 10 


1.03 X 10-' 


0.9022 


0.0942 


2.09 X 10- 


-b 


1.08 xlO- 


b 


4.09 X 10- 


-4 


1.42x10- 


-5 


8.12x10-" 


Dmlt X 3 


1.71 X 10^* 


0.9514 


0.0283 


7.99 X 10- 


-3 


0.0118 




5.13 X 10- 


-5 


8.18x10- 


-8 


1.54x10"** 


Dmlt 


4.60 X lO^*" 


0.9524 


0.0264 


7.73 X 10- 


-3 


0.0134 




4.53 X 10- 


-5 


2.37x10- 


-7 


4.16x10"* 




1.66 X 10^^ 


0.9532 


0.0291 


7.50 X 10- 


-3 


9.71 X 10- 


-3 


3.62 X 10- 


-5 


7.97x10- 


-7 


1.49x10"' 


Dmlt /WO 


4.22 X lO^'' 


0.9506 


0.0328 


7.66 X 10- 


-3 


8.49 X 10- 


-3 


3.18 X 10- 


-5 


1.95x10 


-6 


3.72x10"' 


Dmlt/1000 


4.57 X lO^-* 


0.9515 


0.0327 


6.89 X 10- 


-3 


7.99 X 10- 


-3 


3.06 X 10- 


-5 


2.49x10- 


-5 


3.95x10"* 


Dmlt/10000 


0.0117 


0.9455 


0.0376 


2.53 X 10- 


-3 


1.87 X 10- 


-3 


3.35 X 10- 


-5 


2.27x10- 


-4 


7.42x10"' 



Table 5. Surface abundances (by mass fractions) of the sequences calculated in this work with different mixing efficiencies. All 
sequences correspond to models with a final mass equal to 0. 48545 and Z=0.001. ' This sequence did not undergo the H-flash, 
see text. 



the local flux due to the primary He-flash — and the violent pro- 
ton burning is never sufficient to drive convection further out. As 
a consequence, the "runaway" burning of protons does not de- 
velop and, as can be seen from Table 5, far higher H abundances 
are left than in the standard case. The H-abundance after the pri- 
mary He-flash in this case, may also lead to the conversion of 
He-sdO stars into He-deficient sdB stars, according to Unglaub 
(2005). We note also that a feature of the absence of the runaway 
proton burning (H-flash) is that the convective zone developed 
by the primary HeCF never splits and the surface '^C abundance 
continues to increase until the convective zone merges with the 
outer convection. In this case, the final '^C is much larger than 
the cases in which the H-flash develops. 



4.2. Effect of extra-mixing at convective boundaries 

To analyze the effects of possible extra-mixing episodes, we in- 
cluded extra-mixing zones at each convective boundary follow- 
ing the prescription of Herwig et al. (1997). We therefore al- 
lowed extra-mixing episodes by considering that the mixing ve- 
locities decay exponentially beyond each convective boundary. 
The diffusion coefficient beyond the formal convective bound- 
ary was then given by 

Dov(z) = A)expj-^j (4) 

where Do is the diffusion coefficient of unstable regions close to 
the convective boundary and z is the geometric distance from the 
edge of the convective boundary. We performed simulations with 
three different values of the free parameter / (/=0.005, 0.01 and 
0. 1) around the value inferred from the width of the upper main 
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sequence (/ ~ 0.016) and oxygen abundances of PG1159 stars 
(Herwig 2000). 

Results of these numerical experiments are shown in Table 6. 
In the case of / - 0.005, the final surface abundances are close 
to those of the standard simulations in which no extra-mixing 
was considered. For larger extra-mixing efficiencies H-burning 
becomes more efficient and for values beyond / - 0.01 almost 
no H remains in the star after the H-flash. For values of / < 0.01, 
the evolutionary timescales ('Af, ^At ,^Af; see Table 2 for the 
definition) are unaffected and are similar to the standard case 
(see Table 2). In contrast, for the / = 0.1 case, the time inter- 
val between flashes is significantly reduced (46 d) and the EHB 
evolution is somewhat enlarged by the appearance of breathing 
pulse instabilities at the end of the EHB (212 Myr). The lower 
N abundance in the / = 0. 1 case is because the H is mixed with 
a larger mass of '^C reducing the amount of protons per carbon 
nuclei. 



4.3. Effect of chemical gradients (Vjuj 

Miller Bertolami et al. (2006) recently showed that chemical gra- 
dients CVfj) can affect the development of the rapid ingestion of 
protons into convective zones driven by helium shell flashes. It 
therefore appears important to analyze the possible effect of V// 
in deep-mixing events. To analyze the effect of chemical gradi- 
ents during deep-mixing events, we performed some simulations 
of deep-mixing for a model of 0.4641 Mq and Z=0.02 with the 
theory of convection developed by Grossman et al. (1993), in 
its local version presented by Grossman & Taam (1996^ This 
theory incorporates the effect of V// both in the determination 
of the different mixing regimes and the corresponding mixing 
velocities. 

This set of simulation^ showed that the inclusion of V/i pre- 
vents the complete penetration of the He-core flash convective 
zone into the regions of high H-content. In particular, at the H- 
He transition, the convective zone splits into different convective 
zones and, consequently, the chemical profile at the H-He tran- 
sition acquires a staircase-like aspect. Eventually some of these 
convective zones merge and some H can be burnt. However, the 
main consequence is that in these simulations almost all the H- 
rich envelope is not engulfed and burned. In all of these sim- 
ulations the final abundances remain unchanged after the flash, 
because even the dilution event is prevented due to the yU-baiTier. 

However, the theory adopted here is a local theory and does 
not include any kind of extra-mixing beyond the formal con- 
vective boundaries, such as those produced by overshooting 
(Freytag et al. 1996) or the excitation of gravity waves (Herwig 
et al. 2006, 2007). Our simulations indicate that the inclusion of 
an overshooting efficiency as small as / = 0.001 in our mod- 
els allows the penetration of the yu-barrier and the development 
of the H-flash, the surface abundances that arise in these cases 
are similar to those in which fi effects are not taken into ac- 
count. Lower /-values are difficult to test in our simulations be- 
cause they require very high spatial resolution. Although hydro- 
dynamical simulations of convection in homogeneous mediums 
appears to support the existence of some extra-mixing at convec- 



For simulations adopting the convection theory of Grossman et al. 
(1993), the value of the free parameter a is reduced by V2 from that 
adopted in the standard MLT due to a difference in the definitions of the 
relevant coefficients. 

* We performed several simulations under different numerical as- 
sumptions. 



five boundaries, it should be noted that overshooting is expected 
to be reduced by the //-barrier itself (Canuto 1998). 



4.4. Effect of diffusion in the outer layers 
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Fig. 7. Top: Time evolution of the chemical abundances at r = 
2/3 as a consequence of gravitational settling and radiative lev- 
itation. Initial abundances correspond to one of our DM se- 
quences (Table 2; M= O.48545M0 and Z=0.001) with surface 
parameters Teff - 40000 K and log g - 6.0. Bottom: Chemical 
profile of the outer layers after 10^ yr when the surface was al- 
ready H-dominated. Note the very steep chemical profile close 
to the photosphere (T2/3) which can strongly affect the surface 
values if some mixing process is effective close to the photo- 
sphere. In particular note the supersolar He-abundances below 
T = 50 (rgo), not far from the location where the bottom of the 
convective zone due to the He-ionization should be located. 



We now address the possible effects of the interplay between 
mass-loss, gravitational settling, chemical diffusion, and radia- 
tive levitation in our standard sequences (Tabs. 2 and 3). It is 
well known that the abundances of sdB stars depend strongly 
on the combined effects of diffusion and mass-loss (Unglaub & 
Bues 2001), and the same should be true for He-rich subdwarfs 
(Unglaub 2005). If mass loss is higher than 10"'^Mo/yr, diffu- 
sion would not be able to alter the post-HeCF abundances sig- 
nificantly and the star would remain He-rich during the entire 
EHB evolution. For the low hydrogen abundances of our DM 
sequences, a mass-loss rate of approximately lO^'^^Mo/yr would 
be sufficient to prevent transformation of a He-rich into a H-rich 
star due to the upward diffusion of hydrogen. 
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Table 6. Surface abundances of sequences in which extra-mixing at convective boundaries was considered with different efficiencies 
(/). The sequences correspond to models with a final mass equal to O.4641M0 and Z=0.02. 



From the theoretical predictions of Vink & Cassisi (2002) 
and Unglaub (2008), it follows that for stars with M, k O.SMq, 
T'eff = 40000K, log^ $ 5.5, and solar composition, mass- 
loss rates of M ^ lO'^Mo/yr are expected. The effect of 
diffusion should therefore be negligible. For more compact 
stars with \ogg > 5.5, from the theory of radiatively-driven 
winds, which neglects multi-component effects, weak winds 
with M < lO'^Mo/yr are predicted. However, the results 
of Unglaub (2008) showed that those weak winds cannot be 
"chemically homogeneous" as assumed in previous investiga- 
tions. Hydrogen and helium, at least, can hardly be expelled 
from the star Therefore, with the method described in Unglaub 
& Bues (1998), we performed simulations, in which the absence 
of any mass-loss was assumed. 

In the upper panel of Fig. |7] we show the evolution of the 
surface abundances (at a Rosseland mean optical depth r - 
2/3) predicted by one of our standard DM sequence^ (that 
with 0.48545 and Z - 0.001) as a consequence of gravita- 
tional settling and radiative levitation (at Teg - 40000K and 
logg - 6.0). In the absence of mass loss, the star changes from 
a He-rich subdwarf into a He-deficient one in only ~ 10^ yr. 
Since this timescale is similar to the time it takes the star to 
reach the ZAHB, surface abundances are expected to become 
H-dominated as the star settles onto the ZAHB. This leads to 
the evolutionary channel sketched in Fig. [8], in which He-sdO 
stars become hot sdB stars due to element diffusion. After about 
10^ yr, the surface abundances are close to the values expected 
from the equilibrium condition between gravitational settling 
and radiative levitation obtained by Unglaub & Bues (1998). 
This timescale is significantly shorter than the timescales for 
stellar evolution close to the EHB (65 to 90Myr). For higher 
H/He ratios, as predicted by our shallow-mixing sequences, 
the stellar atmosphere becomes hydrogen-rich on even shorter 
timescales. For the abundances from the SM* sequence with 
O.4725OM0 and Z = 0.02, the equilibrium abundances should 
already be reached by lO^yr. In the absence of mass loss or an- 
other disturbing processes, the surface composition of the ma- 
jority of sdB stars should therefore be given by the equilibrium 
condition between gravitational settling and radiative levitation, 
independently from the evolutionary scenario from which they 
originate. 

It is well known that the equilibrium abundances of helium 
are lower than observed (Michaud et al. 1989). The number ra- 
tios He/H of the hottest sdB stars are approximately 0.01 (e.g. 



' With the exception of the O abundance that was defined to 
be log(A?o/Mot) = -3, instead of log(A'o/A'iot) ~ -5 as follows 
from our simulations, for numerical reasons. For an initial value of 
log(A'o/Mot) ~ ~5, oxygen moves outward in the outer regions and 
its abundance tends to its equilibrium value, which is higher than 

log(A'o/A'iot) 5. In the underlying regions, however, where oxygen 

has helium-like configuration, the outward radiative force is negligibly 

small and oxygen sinks. For an initial abundance of log(A'o/^toi) 5, 

this situation leads to extremely steep concentration gradients in these 
regions where the flow changes its sign and, as a consequence, to nu- 
merical problems. 



Lisker et al. 2005). There is therefore a discrepancy of about 
three orders of magnitude between the predictions and the ob- 
servations. In Fig. |7] we show the predicted abundances as a 
function of depth after 10^ yr We note the very steep chemical 
gradient in the outer regions. Within the stellar atmosphere, the 
helium abundance increases from He/H ^ 10'^ at r = 2/3 to su- 
persolar values at t ^ 50 (see Fig.|7]). We therefore suggest that 
the presence of even weak mixing processes (e.g. due to stellar 
rotation), which tend to smooth concentration gradients, should 
produce higher helium abundances close to the stellar surface. 
Thus the large discrepancy for the helium abundance does not 
necessarily exclude the proposed scenario. In future investiga- 
tions, it should be shown how mixing processes other than mass 
loss would affect the upward diffusion of hydrogen, and the grav- 
itational settling of helium in a H-dominated stellar atmosphere. 
In particular, the possible effect of a thin, outer, convection zone 
due to changes in ionization of helium, predicted by Groth et al. 
(1985) at least for supersolar helium abundances, could affect 
the composition of the stellar atmosphere. For example, in the 
DM sequence analyzed in this section, stellar evolution models 
(without diffusion) show that the bottom of this convective re- 
gion should be located around log(l - m/M*) 12 — with the 

exact location depending on the instantaneous log T^ff and log g 
values — where He abundance is significantly higher (see Fig. 

In the log Tgff - logg diagram of Fig. [8] the line is shown 
above which, according to Unglaub (2008), chemical homoge- 
neous winds may exist for stars of solar metallicity. For stars 
above this line, the mass-loss rates should be sufficiently high 
to prevent diffusion. The chemical composition should then de- 
pend essentially on their evolutionary history. For stars below 
this "wind limit", no chemically homogeneous winds can ex- 
ist. The surface composition of these stars should essentially de- 
pend on atmospheric processes such as gravitational settling and 
radiative levitation. Pure metallic winds may possibly still ex- 
ist, which could produce additional changes in the metal abun- 
dances. Due to the complexity of the situation, no reliable quan- 
titative predictions appear to be presently possible. Our results 
in Fig. |7] indicate, however, that in the absence of mass-loss 
a helium-rich atmosphere may become a hydrogen-rich one, 
even for the low hydrogen abundances predicted by the DM se- 
quences. 

From Fig. [8] it can be seen that our DM sequences cross 
the wind limit before the star settles onto the ZAHB. We then 
expect the onset of diffusion. A similar situation occurs in hot 
white dwarfs on the upper cooling sequence (Unglaub & Bues 
2000). As soon as winds cease, the abundances of heavy ele- 
ments decrease and helium-rich white dwarfs became hydrogen- 
rich. We emphasize that the location of the wind limit in the 
log Teff - logg diagram depends on the chemical composition 
and the line shown in Fig.[8]is for solar abundances. It would be 
necessary to investigate the dependence of the mass-loss rates 
on the abundances of individual elements. From Fig. [8] it can 
be seen that there is at least qualitative agreement between our 
predictions and observational results. The helium-poor sdB stars 
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Fig, 8. Sketch of the proposed channel for the evolution after a DM event as a consequence of chemical diffusion in the absence 
of winds (thick arrows mark the sense of evolution). A similar situation would be expected for post-SM/SM* sequences with weak 
winds (M ~ lO^'^Mo/yr). Upper panel: proposed evolution in the abundances-Teff plane, abundances taken from Stroer et al. 
(2007) and Lisker et al (2005). Lower Panel: Evolution in the g - Teff plane. Continuous lines correspond to the evolution as a 
He-rich subdwarf, while dotted lines mark the evolution after the ZAHB when diffusion would have turned the atmosphere into a 
H-dominated one. Black dashed lines show the wind limit for solar-like abundances and for two different metallicities (thick line: 
Zq, thin line: Zq/3). [color figure only available in the electronic version] 



are preferably in the region below the wind limit, where grav- 
itational settling could be effective. The He-rich subdwarfs are 
preferably above the wind limit, where we expect no diffusion. 
The composition of these stars should show clear signatures of 
their evolutionary history. 

5. Discussion 

In Sect. 3.3, we showed that standard sequences lead to final 
H abundances (see Table 3) immediately after a DM event that 
are much lower than inferred for He-sdO stars by Stroer et al. 
(2007). Only a minority of standard sequences (those which un- 
dergo SM/SM* events) display abundances similar to He-sdO 
stars immediately after the flash. Also, as we have already men- 
tioned, the predicted location in the Teff would differ from that 
observed. In particular, since the EHB evolution is more than 40 
times longer than the evolution from the primary HeCF to the 
ZAHB, most of the stars should be clustered around the EHB, 
which is not observed. However, there are some viable scenar- 
ios in which our results are consistent with observations. These 



scenarios depend on the existence or not of weak homogeneous 
winds in He-rich subdwarf stars as we now discuss. 

As mentioned in Sect. 4.4, in a very recent study Unglaub 
(2008) showed that weak homogeneous winds are impossible in 
hot subdwarfs at high gravities. In this case, as shown in Fig.|7] 
the low amount of H remaining in our standard sequences after a 
DM event is sufficient to convert a He-rich to a He-deficient sur- 
face composition in about 10^ yr Within this picture, after the 
DM event the star first evolves towards the ZAHB star as a He- 
rich subdwarf. Then, when gravity has increased sufficiently, the 
star crosses the wind limit causing homogeneous winds to halt, 
then diffusion leads to an increase in the H abundance as the star 
evolves towards the ZAHB (both processes occur on timescales 
of the order of 1 Myr). We note in particular (Fig.|7]) that in only 
a few thousand years the H abundance rises from its post-DM 
value to the values usually measured in He-sdO stars (Stroer 
et al. 2007) and thus the correct H-abundances are predicted. 
When the HB has been reached the evolutionary timescale be- 
comes more than a factor of 40 longer (65 to 90 Myr, see tab. 
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2) and the star turns into a H-dominated sdB star This possible 
evolutionary channel is sketched in Fig. [8] In this case no clus- 
tering of He-sdO stars should be expected close to the EHB, in 
good agreement with the results of Stroer et al. (2007). It is inter- 
esting to note that the effective temperatures expected for these 
stars, when on the HB (and with a H-dominated atmosphere), 
agree with those of the hotter sdB stars (which also have higher 
He abundances than their cooler counterparts). The conversion 
of surface abundances with C>N into surfaces with N>C as H 
diffuses outwards, might also be related (see Figs. |7] and [8]l to 
the lack of C-rich He-sdO stars with lower Nue/Nu values in the 
study of Stroer et al. (2007). Within this scenario and in the ab- 
sence of any other channel that creates hot-subdwarfs, the num- 
ber ratio of He-sdO stars to hot (Teff > 33 500) sdB stars should 
be of the order of 1/40 in a complete sample (although higher 
ratios can still be possible if some process that disturbs diffu- 
sion is present). Once the star finishes the He-core burning stage 
(or EHB), it moves above the wind limit again (see Fig. [8]) and 
homogeneous winds again become possible. If this is so, since 
the H-rich layer is very thin, a new conversion into a He-sdO star 
might happen (now at higher surface gravities and temperatures). 
Finally, it is interesting to note in Fig. [T] that the C-abundance 
decrease rapidly by about a factor of 25 and the star changes its 
surface composition from C>N to C<N. This occurs while the 
star is still strongly H-deficient. We note that, at this point (be- 
tween t=1000 and 10000 yr in Fig.|7]), the star displays surface 
abundances qualitatively similar (H~ 2.5 x 10"^, N~ 2 x lO^^* 
and C<N by mass fraction), although quantitatively different, to 
those inferred by Lanz et al. (2004) for the peculiar star LB 1766. 

In the case of post-SM/SM* stars our simulation indicates 
that the atmosphere should turn into a H-dominated one in about 
1000 yr from the moment in which homogeneous winds cease. 
Then, if homogeneous winds do not exist at high gravities, post- 
SM/SM* stars would only be He-rich while they are above the 
wind limit (see Fig. [8]). This is clearly not the case for most He- 
sdO stars in the sample of Stroer et al. (2007) but it does describe 
the observed surface properties of JL87 — reff = 29 000 K and 
log^ = 5.5, Lanz et al. (2004); T^ff = 26 000 K and logg = 4.8, 
Ahmad et al. (2007). 

We emphasize that the possible conversion of He-sdO stars 
into hot sdB stars in the HB discussed in the previous paragraphs 
is interesting, not only because it predicts the correct distribution 
in the log g - Teff diagram but also in view of the significantly 
higher He abundances of sdB stars with T^g > 33 500K com- 
pared with their cooler counterparts, which might reflect differ- 
ences in their previous evolution. 

In contrast, if the effects of element diffusion were of minor 
importance for all He-rich subdwarfs, then only the abundances 
resulting from SM/SM* events would be consistent with the sur- 
face properties measured for He-sdO stars (Stroer et al. 2007), 
unless DM abundances differ strongly from those predicted by 
the standard MLT theory (see Sect. 4.1). We note that the frac- 
tion of stars undergoing SM/SM* events corresponds to only a 
small fraction of hot-flashers in standard sequences. In this con- 
text, our results in Sect. 4.3 indicate that the role of chemical 
gradients should be analyzed to determine whether the fraction 
of shallow-mixing episodes might be increased by the effect of 
chemical gradients that tend to prevent the penetration of con- 
vection into the H-rich envelope. 

The hot-flasher scenario therefore can qualitatively predict 
the correct surface properties and distribution of He-sdO stars. It 
is then interesting to note that the lack of close binaries (with pe- 
riods of the order of days as observed for sdB stars) in He-sdO 
stars is also expected naturally within the hot-flasher scenario. 



Very close binaries with final periods of the order of days will 
indeed fill their Roche lobes before reaching the RGB-tip and 
never undergo the HeCF. Also, it is not clear if it would be pos- 
sible to have hot-flashers within a common envelope system. In 
any case, in order to have a hot-flasher within a close binary sys- 
tem, a very fine tuning of the initial masses and periods seems 
to be needed. Consequently, post hot-flasher H-deficient stars 
should not be usual in these systems. The lack of close binary 
systems with He-sdO components should then not be taken as 
an argument against the hot-flasher and in favor of the merger 
scenario. In contrast, wide systems would be a much better way 
of distinguishing between both possible evolutionary origins for 
He-sdO stars, because they should be more probable within the 
hot-flasher picture. 

6. Final Remarks 

We have performed ID stellar evolution simulations of the hot- 
flasher scenario for a wide variety of cases and metallicities. 
Based on these results, we have given a detailed description of 
the surface characteristics of standard sequences. For these se- 
quences, we have studied how abundances could be altered by 
element diffusion processes. We have then extended the scope 
of our work by studying deviations from standard assumptions. 
Our results can be summarized as follows: 

- We have carried out an ample exploration of the parameter 
space of the hot-flasher scenario that will allow far more ef- 
fective tests of its predictions. Our results confirm the partial 
results by Cassisi et al. (2003) and Lanz et al. (2004) for the 
deep-mixing and shallow-mixing events respectively. 

- Our deep-mixing sequences have abundances similar to 
those recently inferred for PG1544+488, confirming that this 
star can be the result of a DM event. Our results also con- 
firm the result of Lanz et al. (2004) that JL87 abundances are 
compatible with a recent SM event for this star None of our 
sequences show abundances similar to those inferred for LB 
1766 and this object is therefore probably not the result of a 
recent hot-flasher event. However, it is worth noting that af- 
ter a DM event, element diffusion might produce abundances 
qualitatively similar to those of LB 1766. 

- Standard hot-flasher sequences predict C/N ratios in the 
range from ~ 9.5 to ~ 0.85. This is particularly interesting in 
the light of the N-rich and C-rich classes described by Stroer 
et al. (2007), although no clear correlation is found between 
those classes and our sequences in the g-T^s diagram. In the 
absence of weak winds, however, the lack of C-rich He-sdO 
stars with lower He/H values in the sample of Stroer et al. 
(2007) can be qualitatively understood in terms of element 
diffusion. 

- The location of our sequences in the g - Teff diagram show 
qualitative agreement with the parameters inferred by Stroer 
et al. (2007) for He-sdO, while they predict (as already noted 
by Stroer et al. 2007) that these stars should cluster around 
the EHB, which is not observed. An interesting (and plau- 
sible, in view of the results of Unglaub 2008) possibility for 
this discrepancy to be removed is if He-sdO stars evolve from 
a region in the g - T^ff diagram where homogeneous winds 
can exist, into a region where such winds are impossible. 
This should lead to the onset of H-diffusion which, as we 
have showed in Sect 4.4, allows the conversion of He-sdO 
stars close to the ZAHB into H-dominated sdB stars even for 
the low abundances resulting from DM events. This scenario 
is extensively discussed in Sect. 5. 
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In the presence of weak homogeneous winds only the surface 
properties arising from SM/SM* episodes would be consis- 
tent with those inferred for He-sdO stars, unless DM abun- 
dances differ strongly from those predicted by the standard 
MLT with diffusive convective mixing. 

- In all cases in which a conversion of He-sdO into hot H- 
rich sdB stars is possible, we have shown that the hot-flasher 
scenario reproduces correctly the observed properties (abun- 
dances and their distribution in the T^ff - g plane) of He-sdO 
stars and links them, as immediate progenitors, to the hot- 
ter sdB stars. This is interesting in view of the higher He- 
abundances of the hotter sdB stars. 

- We have argued that the lack of close binaries among He-sdO 
stars cannot be used as an argument against the hot-flasher 
scenario because this dearth is also naturally expected within 
this scenario. 

- Our study of departures from standard assumptions shows 
that the role of chemical gradients as an extra barrier to the 
penetration of the inner convective zone into the H-rich en- 
velope has to be studied in detail. This is particularly im- 
portant in determining if the fraction of SM/SM* episodes 
is higher than predicted by standard sequences. In contrast, 
extra-mixing processes alone does not seem to produce sig- 
nificant departures from standard predictions. 

We conclude that the hot-flasher scenario is a viable expla- 
nation of the formation of He-rich subdwarfs and that the under- 
standing of some key points needs to be clarified before a final 
statement can be made. Special attention should be paid to ana- 
lyzing the existence and intensity of disturbing processes, such 
as winds and turbulence, in the outer layers of He-rich subd- 
warfs. Also, the search for wide binaries among He-sdO stars is 
necessary to distinguish between the predictions of the merger 
and hot-flasher scenarios. Finally, a better understanding of the 
burning and mixing of H in the hot interior (probably by means 
of hydrodynamical simulations of the event) is badly needed. 
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Note Added in Proof: After acceptance, the authors became 
aware of the close binary nature of PG 1544+488 (Ahmad, A. 
& Jeff'ery C. S. 2007, ASPC, 391, 261). This fact raises serious 
doubts on a possible deep-mixing origin for this star 
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